Angular magnetoresistance oscillations in quasi-one-dimensional organic conductors in 

the presence of a crystal superstructure 
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We study the effect of crystal superstructures produced by orientational ordering of the Re04 and 
CIO4 anions in the quasi-one-dimensional organic conductors, (TMTSF) 2 Re0 4 and (TMTSF) 2 00 4 , 
on the angular magnetoresistance oscillations (AMRO) observed in these materials. Folding of the 
Brillouin zone due to anion ordering generates effective tunneling amplitudes between distant chains. 
These amplitudes cause multiple peaks in interlayer conductivity for the magnetic-field orientations 
along the rational crystallographic directions (the Lebed magic angles). Different wave vectors of 
the anion ordering in (TMTSF) 2 Re0 4 and (TMTSF) 2 C10 4 result in the odd and even Lebed an- 
gles, as observed experimentally. When a strong magnetic field is applied parallel to the layers and 
perpendicular to the chains and exceeds a certain threshold, the interlayer tunneling between differ- 
ent branches of the folded electron spectrum becomes possible, and interlayer conductivity should 
increase sharply. This effect can be utilized to probe the anion ordering gaps in (TMTSF) 2 C10 4 
and (TMTSF) 2 Re0 4 . An application of this effect to k-(ET) 2 Cu(NCS) 2 is also briefly discussed. 

PACS numbers: 74.70.Kn, 72.15.Gd, 73.21.Ac 



I. INTRODUCTION 

The quasi-one-dimensional (Q1D) organic conductors 
(TMTSF) 2 X (where TMTSF is tetramethyltetraselena- 
fulvalene and X represents a monovalent anion, such 
as PF6, CIO4, or Re04) have very interesting physical 
properties, including the quantum Hall effect and pos- 
sibly triplet superconductivity 1,2 . These materials con- 
sist of parallel conducting chains along the x axis, ar- 
ranged in layers with the interchain spacing b along the 
y axis and the interlayer spacing c along the z axis. 
The electron-tunneling amplitudes between the TMTSF 
molecular sites are highly anisotropic in the three direc- 
tions: t a :t b :t c = 2500 : 250 : 10 K 1 . 

These materials exhibit the angular magnetoresis- 
tance oscillations (AMRO), where resistivity strongly 
changes as a function of the magnetic-field orientation. 
There are three basic types of AMRO: the Lebed magic 
angles 3-7 for the magnetic-field rotation in the (y, z) 
plane, the Danner-Kang-Chaikin (DKC) oscillations in 
the (x,z) plane 8 ' 9 , and the third angular effect in the 
(x,y) plane 10-12 . The Lebed oscillations manifest them- 
selves as sharp peaks in the interlayer conductivity a zz 
occurring when the magnetic field points from one chain 
to another along a rational crystallographic direction, as 
illustrated in Fig. 1. Approximating the triclinic crystal 
lattice of (TMTSF^X by the orthogonal one, the magic 
Lebed angles can be written as 



fects coexist. They found that the Lebed oscillations are 
enhanced when B x ^ 13 , and the DKC oscillations still 
exist in the presence of B y 7^ 14 . 

Although initially the different types of AMRO were 
treated as separate phenomena, a unified picture emerged 
in the recent years due to substantial experimental and 
theoretical progress. A three-dimensional visualization 
of the experimentally measured a zz (B) (Ref. 15) demon- 
strated that the different types on AMRO can be viewed 
as modulations of the basic Lebed resonances. Measure- 
ments with carefully placed electric contacts 16 proved 
that AMRO exist only in the transverse resistance R zz 
and not in the longitudinal resistance R xx along the 
chains. Theory always predicted this difference, but 
many experiments observed AMRO in R xx as well be- 
cause of the mixing between different components of the 
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where n and m are integer numbers, and B = 
(B x , B y ,B z ) = B (sin 9 cos if, sin 9 sin tp, cos 9) is the mag- 
netic field. Experimentally, the Lebed effect is the most 
pronounced for m = 1. Lee and Naughton 13,14 studied 
AMRO for generic orientations of B, where all three ef- 



FIG. 1: (Color online) A view along the chains of a Q ID metal 
with the anion ordering at a wave vector Q. The filled and 
open circles represent the chains with the energies ±E g . (a) 



(TMTSF) 2 Re0 4 and Q 
and Q = (0,1/2,0). 



(0,1/2,1/2). (b) (TMTSF) 2 C10 4 
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conductivity tensor. AMRO were found not only in the 
dc conductivity, but also in the ac conductivity at mi- 
crowave frequencies 17 ' 18 . The ac measurements were in- 
terpreted in terms of the so-called period orbit resonance 
(POR) 19 , which is a generalization of the cyclotron reso- 
nance to more complicated (e.g., open) Fermi surfaces 20 . 
The ac resonances occur at the angles depending on fre- 
quency w and deviating from Eq. (I) 17 ' 18 , so the Lebed 
magic angles are not truly magic 21 ' 22 . This observa- 
tion eliminates theoretical scenarios proposing a radical 
change in the ground state of the system depending on 
the magnetic-field orientation along the magic or non- 
magic angles. This conclusion is also supported by the 
absence of any angular effect in NMR 23 . 

Given these experimental facts, AMRO most likely 
represent some sort of a resonance effect in the dc and ac 
transport coefficients. The first theoretical calculation 
along these lines was done in Ref. 24 using the Kubo 
formula with the electron wave functions for a magnetic 
field in the (y, z) plane. This quantum-mechanical calcu- 
lation was then generalized to include the B x component 
of the magnetic field 25 and the anion superstructure of 
(TMTSF) 2 C10 4 26 ' 27 . In another theoretical approach, 
the Boltzmann kinetic equation was solved for a con- 
stant relaxation time r by using quasiclassical electron 
trajectories on the Fermi surface 11-14 ' 19,20 ' 28-31 . This so- 
lution can be written in a general form using the so-called 
Shockley tube integral 32 or the Chambers formula 33 , see 
also the book 34 . In the third theoretical approach, the 
interlayer conductivity was calculated using a perturba- 
tion theory in the electron-tunneling amplitude between 
two layers 35-37 . In this approach, AMRO originate from 
the Aharonov-Bohm quantum interference in interlayer 
tunneling in the presence of a magnetic field 37 . All these 
three seemingly different theoretical approaches produce 
the same final results and are essentially equivalent. 

Despite substantial progress in understanding of 
AMRO in Q1D conductors, some experimental results 
remain unexplained. One open problem is the angu- 
lar oscillations of the Nernst effect 38 . Another unre- 
solved problem is the angular minimum and saturation 
of the interlayer resistivity R zz observed for a magnetic 
field in the y direction 6 ' 14 ' 15 ' 39 . Although the manifes- 
tations of AMRO arc qualitatively similar in all mem- 
bers of the (TMTSF)2X family, direct comparison of the 
measurements in (TMTSF) 2 PF 6 , (TMTSF) 2 C10 4 , and 
(TMTSF) 2 Re0 4 shows substantial differences 39 . 

For a magnetic-field rotation in the (y, z) plane with 
B x = 0, only three strong Lebed peaks in a zz with n = 
0, ±1 are observed in (TMTSF) 2 PF 6 6 - 39 . When special 
care is taken to ensure that B x = 0, the very weak peaks 
with n = ±2 in (TMTSF) 2 PF 6 disappear completely 15 . 
In contrast, in (TMTSF) 2 Re04, strong Lebed oscillations 
are observed up to n = ill 40 . In (TMTSF) 2 C10 4 , the 
Lebed oscillations are much weaker in amplitude than in 
(TMTSF) 2 PF 6 and (TMTSF) 2 Re0 4 39 , but many Lebed 
resonance can be detected after differentiation of the data 
with respect to the angle of rotation 4,5 . The strength 



of the DKC oscillations is also very different in these 
materials. The DKC oscillations are quite strong in 
(TMTSF) 2 C10 4 , where they were originally discovered 8 . 
In (TMTSF) 2 PF 6 , Ref. 9 found very weak DKC oscilla- 
tions, but Ref. 39 found them to be substantial. However, 
in (TMTSF) 2 Re0 4 , the DKC oscillations are extremely 
weak and almost invisible 39 . This dramatic difference in 
manifestations of AMRO in the three materials requires 
a theoretical explanation. 

When a magnetic field is rotated in the (y, z) plane 
at B x = 0, the theoretical calculations cited above show 
that the Lebed peaks in a zz can exist only for those magic 
angles (n, m) where the interchain tunneling amplitudes 
in the directions nb + mc are present 24,41 . It is reason- 
able to expect that the interplane tunneling amplitudes 
in (TMTSF) 2 PF 6 exist between the nearest and next- 
nearest chains in the c and c ± b directions (see Fig. 
1). This would explain why only the Lebed resonance 
with n = 0, ±1 are observed in (TMTSF) 2 PF 6 . How- 
ever, many magic angles with big numbers n are observed 
in (TMTSF) 2 C10 4 and (TMTSF) 2 Re0 4 . It is hard to 
imagine that direct electron overlap exists between the 
chains separated by 11 interchain distances. 

One way to resolve this problem is to take into ac- 
count the nonlinear electron dispersion along the chains. 
(All theoretical papers cited above make a linearized ap- 
proximation for the electron dispersion along the chains.) 
The first attempt in this direction was made in Ref. 42, 
and a more systematic study was presented in Refs. 43 
and 44. The nonlinearity can indeed generate an effect 
similar, albeit not completely equivalent, to the presence 
of many interchain tunneling amplitudes. However, the 
nonlinearity alone is not sufficient to explain the dif- 
ferences in AMRO between the three compounds. An- 
other problem is the absence of the DKC oscillations in 
(TMTSF) 2 Re0 4 . One might think that quantum co- 
herence is too low in this material, but the existence 
of 21 Lebed oscillations clearly refutes this idea 40 . We 
see that a detailed theoretical understanding of AMRO 
in the (TMTSF) 2 X materials is challenging and requires 
additional ideas. 

We believe that the key to understanding the differ- 
ences in AMRO is the presence of anion ordering in 
(TMTSF) 2 C10 4 and (TMTSF) 2 Rc0 4 and its absence in 
(TMTSF) 2 PF 6 . PF 6 is an octagonal centrosymmetric 
anion, which does not experience any orientational order- 
ing at low temperatures. In contrast, C10 4 and Re0 4 are 
tetragonal anions without inversion symmetry. Because 
their crystal sites have inversion symmetry, these anions 
have two different orientations of the same energy. At low 
temperatures, the anions experience orientational order- 
ing and produce crystal superstructures 1 with the wave 
vectors Q = (0,1/2,0) in (TMTSF) 2 C10 4 (under ambi- 
ent pressure) and Q = (0,1/2,1/2) in (TMTSF) 2 Re0 4 
(under pressure greater than about 10 kbar), as shown 
in Fig. 1. Formation of a crystal superstructure affects 
electron spectrum by folding the Brillouin zone. In this 
paper, we show that reconstruction of the electron dis- 
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persion caused by the anion ordering generates effec- 
tive tunneling amplitudes between many distant chains. 
This effect explains why many Lebed angles are ob- 
served in (TMTSF) 2 Re0 4 and (TMTSF) 2 C10 4 , but not 
in (TMTSF) 2 PF@. ft also explains why the magic angles 
[Eq. (1)] are observed only for odd n in (TMTSF) 2 Re0 4 
(Rcf. 40) and only for even n in (TMTSF) 2 C10 4 (Refs. 4 
and 5) at m = 1. We also explain the differences in the 
DKC oscillations within the same framework. 

In contrast to the previous theories of AMRO for 
the anion superstructure of (TMTSF) 2 C10 4 26 ' 27 ' 30 ' 43 , we 
take into account the direct effect of anion ordering on 
the interlayer tunneling amplitude, which is especially 
important for (TMTSF) 2 Rc0 4 . In this way, we can cap- 
ture the characteristic features of AMRO in the three 
compounds without invoking the nonlinearity of the lon- 
gitudinal electron dispersion 43,44 . 

In the second part of the paper (Sec. VI), we study 
the effect of a strong magnetic field parallel to the layers. 
We show that, when B y is strong enough and exceeds a 
certain threshold related to the anion gap E g , the inter- 
layer tunneling between different branches of the folded 
electron dispersion becomes possible, and a zz should in- 
creases sharply. Experimental observation of this effect 
would allow direct measurement of E g . This effect can 
be also applied to study the interband tunneling in n- 
(ET) 2 Cu(NCS) 2 . A theory of this effect cannot be formu- 
lated within the framework of quasiclassical orbits on a 
warped Fermi surface. We calculate an interlayer conduc- 
tivity in the presence of anion ordering using the quan- 
tum limit, where the electron wave functions are confined 
to the layers due to a strong parallel magnetic field 45-48 . 



II. CALCULATION OF INTERLAYER 
CONDUCTIVITY 

The general form of the electron dispersion in a Q1D 
metal is 



e(fc) = ±hv F (k x =F k F ) +e±{k y ,k z ), 



(2) 



where the energy e is measured from the Fermi energy, 
and k — — {kx i &y i 

k z ) is the electron wave vector. Here we 
linearize the dispersion along the chains with the Fermi 
velocity vf near the Fermi wave vectors ±k F . There 
are two sheets of the open Fermi surface, but we present 
calculations only for the sheet with +vp. Since ( c « (i, 
we can expand the transverse dispersion e± to the lowest 
order in the interlayer tunneling amplitude t c , 

e±(k y ,k z ) = 2t b £ v {k v b) + 2t c f(k y b) cos(k z c). (3) 

For a simple model with electron tunneling between 
the nearest chains in the absence of a superstructure, 
Eq. (3) reduces to a standard tight-binding expression 
with Sy(ky) — cos(kyb) and f(k y b) = 1. However, we 
will show in Sees. IV and V that a nontrivial function 
f(k y b) appears in the interlayer tunneling term in the 



presence of anion ordering. This effect was not consid- 
ered in previous literature and plays a crucial role in our 
consideration. 

From the dispersion relation (2), we obtain the electron 
velocity v = de/hdk, 

v x =v F , v v~-^-[jr-> v * = —f{k y b)sm{k z c). 



(4) 

In the quasiclassical approximation, the time-dependent 
electron wave vector feW follows the equation of motion, 



h—— = ev {t) x B, 
dt 



(5) 



where e is the electron charge, and the magnetic field B 
is in the SI units. Given that v x — v F » v z , we find 



dky^ evpB z 



fc(*>--^ + fc(°) oj - ebVpBz (6) 

where u> c is the analog of the cyclotron frequency for the 
open Fermi surface. The equation of motion for k z is 



dk z = ^ ( v F Bydt - 2tb ^ x ^£lL di 



h dk v 



(7) 



Using dky/dt from Eq. (6), we get 

ck^^B^t + B'^yik^+ck^, 
where we introduced the dimcnsionless parameters 



„/ = ^c , _ B x 2t b c 
y B z b' x B z Uv F 



(9) 



The variables B' y and B' x are proportional to the tangents 
of the magnetic field projections onto the (y, z) and (x, z) 
planes, respectively. 

The interlayer conductivity a zz is given by the Shock- 
ley tube integral 34 , 



4e 2 ffdk^dk^ 



h JJ (2n) 3 v F 



dtv z {k^)v z (k^)e t{1 



Jr — iuj) 



(io) 

where r is a relaxation time, and the factor 4 comes from 
the two spin projections and the two sheets of the Fermi 
surface. Substituting Eqs. (4), (6), and (8) into Eq. (10), 
we find the real part of o zz 

2-rr oo 

tt z Ti u) c VFb rr J J 
+ o o 

x cxp{iB' x [e v ((t>) - e v {<j) + r])} - r?[l/w c r - iB' y =F iw/w c ]}, 

where <j) = bk y and r\ = —uj c t. Expanding the periodic 
functions /(</>) e lB * £ ^ in Eq. (11) into the Fourier series 
with the coefficients 

MK) = ^J e-^me^vW # , (i2) 
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we obtain 



\MB' X ) 



+ (w c T) 2 (n — B' =F lu/uj c ) 2 



(13) 



Here er = (4e 2 i 2 -rc)/(7r?i 3 v,F6) is the interlayer dc con- 
ductivity at B = 0, and the ± terms are the contributions 
from the two sheets of the Fermi surface. In the rest of 
the paper, we shall focus on the dc conductivity a zz at 
uj = 0, although Eq. (13) also gives the ac conductivity. 

The Lebed effect corresponds to the resonant peaks of 
a zz in Eq. (13) achieved at B' = n, where the condition 
(1) for to = 1 is satisfied. In a simple model without 
anion ordering, where e y — cos(k y b) and / = 1, Eq. 
(12) reduces to A n (B' x ) = i n J n (B' x ), where J„ is the 
Bessel function. In this case, Eq. (13) reproduces the 
result found in Refs. 25,31,35-37. However, the coeffi- 
cients J n (B' x ) vanish for n ^ at B x = 0, so there are 
no Lebed oscillations in this model for a magnetic-field 
rotation in the (y, z) plane. The DKC effect originates 
from the oscillations of J n (B' x ) vs B' x in the numerator 
ofEq. (13). 

Interestingly, Eq. (13) with \A n \ 2 = J%{B' X ) and to = 
is exactly the same as the equation 49 ' 50 that describes the 
Mach-Zehnder interference in a superconducting qubit 
driven by an ac electric field and subjected to a dc 
bias 49-53 . The two states of the qubit correspond to the 
two adjacent layers of a Q1D conductor coupled by the 
tunneling amplitude t c . The frequency of the ac field for 
the qubit maps to the frequency uj c in Eq. (6), the de- 
tuning of the qubit maps to B' y td c — ecvpBy/h, and the 
amplitude of the ac modulation maps to B' x in Eq. (9). 
The contour plot of Eq. (13) shown in Fig. 2 of Ref. 37 
is exactly the same as in Refs. 49 and 50, and it repre- 
sents the so-called Bessel staircase. The same equation 
also appears in the theory of laser cooling in ion traps . 
This correspondence is not just a mathematical curios- 
ity, but it also reflects profound similarity between these 
highly coherent quantum system, where the oscillatory 
patterns are caused by phase interference due to applied 
electric and magnetic fields. 



III. INTERLAYER CONDUCTIVITY IN 
(TMTSF) 2 PF 6 WITHOUT ANION ORDERING 



transverse dispersion relation can be written as 

s±(k y , k z ) = 2ti, cos(kyb) + 2 ti cos(fc z c + lk y b). (15) 

i 

Equation (14) is the special case of Eq. (15) with t = t c 
and t±i = t' c . 

Generalizing the derivation presented in Sec. II to the 
transverse dispersion relation (15), we find that the in- 
terlayer conductivity a zz is given by Eq. (13) with the 
following coefficients A n (Ref. 55) 



A n {B' x ) = ^Y, in+lt l J n+l{K)- 



(16) 



In the case of (TMTSF) 2 PF 6 , Eqs. (12) and (14) or 
Eq. (16) give 



A n (B' x ) = %^J n (B' x )+i^J^ 1 (B' x )+i n - l "fJ n -i{B' x )- 

(17) 

Substituting Eq. (17) into Eq. (13), we obtain o zz for 
(TMTSF) 2 PF 6 . When B' x = 0, Eq. (17) gives non-zero 
coefficients A n only for n — and n = ±1. Thus, Eq. (13) 
exhibits the Lebed peaks only at n = and n = ±1 with 
the heights proportional to t\ and (t' c ) 2 for a magnetic- 
field rotation in the (y, z) plane. 

When wc consider the DKC oscillations at B' y = 0, 
i.e., for a magnetic-field rotation in the (x, z) plane, the 
sum in Eq. (13) is dominated by the term with n = 0, 
because the other terms have the big factor (w c r) 2 in the 
denominator. Keeping only the term with n = and 
using Eq. (17), we can write approximately, 



-A 



TO 



J (B' x ) + 2ifj 1 (B' x ) 



(18) 



When t' c = 0, Eq. (18) vanishes for the angles where 
Jo(B' x ) = 0, which is a manifestation of the DKC oscilla- 
tions. However, in the presence of t' c ^ 0, Eq. (18) docs 
not vanish for any angles, so the DKC oscillations are 
partially suppressed, although some modulation of a zz 
vs. B' x remains. We see that the presence of tunneling 
amplitudes t\ to more distant chains enhances the Lebed 
oscillations but suppresses the DKC oscillations. This 
conclusion was already made in Ref. 37. 



Let us first discuss the case of (TMTSF) 2 PF 6 , which 
does not have anion ordering. In order to observe more 
than one Lebed angle, we need to introduce the tunneling 
amplitude t' c between next-nearest neighboring chains, as 
shown in Fig. 1(b). Including this term in the transverse 
dispersion (3), wc find for (TMTSF) 2 PF 6 

t> 

s y (<j>) = cos</>, f{4>) = 1 + 2— cos 4>, 4> = bk y . (14) 

In a more general case, where the amplitudes t n corre- 
sponding to the tunneling vectors c + nb are present, the 



IV. ANION ORDERING IN (TMTSF) 2 Re0 4 

The Re04 anions order with the wave vector Q = 
(0, 1/2, 1/2) under pressure. This causes the energies of 
the odd and even chains to split by ±E g , as illustrated 
in Fig. 1(a). The Hamiltonian of interchain tunneling is 
described by a 2 x 2 matrix representing the even and 
odd chains 56 : 



Eg 2tb cos(kyb) + 2t c cos(fc z c) 



'9 

c.c. 



-E n 



(19) 



5 



1 (a) 
0.5 


-0.5 



-1 




-1.5 



-0.5 

(|>/27C 



0.5 



1.5 



FIG. 2: (Color online) (a) Plot of the function f((f>) given 
by Eq. (22). (b) Plot of the second term in f(<f>) given by 
Eq. (25). In both plots, E g /2t b = 0.1. 



The eigenvalues of the matrix (19) give the transverse 
electron dispersion relation, 



e_L = ±J[2t b cos(k y b) + 2t c coa(k z c)} 2 + E 2 . (20) 

Expanding Eq. (20) to the zeroth and first order in t c , 
we find the functions e y (k v ) and f(k y ) in Eq. (3) 



e y (<j>) = ±^ 'cos 2 </> + (E g /2t b y, 
cos 6 



m = ± 



^cos 2 <P + (E g /2t b y 



bk y , (21) 
(22) 



The function /(</>) in Eq. (22) is close to a square wave 
for Eg/tf, <C 1, as shown in Fig. 2(a). Its Fourier coeffi- 
cients A n , given by Eq. (12) with B' x = 0, are non-zero 
only for odd n and decay as 1/n. Transforming Eq. (3) 
from the momentum space to the real space, we find that 
the Fourier coefficients of f(k y b) generate effective in- 
terplane tunneling amplitudes along the vectors c + nb 
with odd n, which are shown in Fig. 1(a) by the arrows. 
Initially, the model has only the tunneling amplitudes 
t b and t c between the nearest chains, but the anion or- 
dering generates effective tunneling amplitudes between 
many chains. The higher-order expansion of Eq. (20) in t c 
would generate effective tunneling amplitudes along the 
vectors mc+nb with m and n of the same parity between 
the sites of the same type, either open circles or closed 
circles in Fig. 1(a). However, one should keep in mind 
that this heuristic real-space picture is an oversimpli- 
fication, and an accurate calculation in the momentum 
space should be performed as described above. 

In Fig. 3 we show the normalized dc conductivity cal- 
culated from Eq. (13) for B' x = and uj c t = V50 us- 
ing the Fourier coefficients A n from Eq. (12). Since 
A n only for odd n, therefore u zz has peaks only at 
the odd Lebed angles, as shown in Fig. 3 and observed 
in (TMTSF) 2 Re0 4 4 °. The higher-order expansion of Eq. 
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FIG. 3: (Color online) Normalized interlayer conductivity 
a zz /a calculated from Eq. (13) for (TMTSF) 2 Re04 and plot- 
ted vs. By at B' x — 0, shown in the linear (left) and logarith- 
mic (right) scales. 



(20) in t c would generate peaks at the Lebed magic angles 
with m and n of the same parity in Eq. (1), as observed 
in Ref. 40. Because of the anion superstructure, Eq. (20) 
is highly non-linear in cos </>; so its Fourier expansion gen- 
erates a big number of harmonics, which produce a big 
number of Lebed peaks in AMRO. This is the qualita- 
tive reason why so many Lebed peaks are observed in 
(TMTSF) 2 Re0 4 , in contrast to (TMTSF) 2 PF 6 , which 
has no anion superstructure. 

Figure 4 shows a contour plot of ln(a zz /ao) vs. B' x and 
B' y , as calculated from Eq. (13) using Eqs. (12), (21), 
and (22). The conductivity is maximal at the vertical 
lines corresponding to the odd Lebed magic angles. At 
a fixed Lebed angle, the weak modulation of a zz vs. B' x 
(along a vertical line) corresponds to the DKC oscilla- 
tions. Figure 4 shows that the DKC oscillations are very 
weak, because the coefficients A n (B' x ) [Eq. (12)] do not 
have zeros vs. B' x in the presence of anion ordering, un- 
like the Bessel functions J n (B' x ) in a simple model. This 
is a theoretical explanation of why the DKC oscillations 
in (TMTSF) 2 Re04 are very weak and barely detectable 
experimentally 40 . 



V. ANION ORDERING IN (TMTSF) 2 C10 4 

In the case of (TMTSF) 2 C104, in order to observe mul- 
tiple Lebed angles, we need to take into account the tun- 
neling amplitude t' c introduced in Sec. Ill and shown in 
Fig. 1(b). For the anion ordering with Q — (0,1/2,0), 
the interchain tunneling is described by the Hamiltonian, 



Eg + 2t c cos(fc z c) cos(kyb) [2t b + it' c cos(fc z c)] 
c.c. — Eg + 2t c cos(k z c) 



(23) 
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FIG. 4: (Color online) Contour plot of ln(cr zz /<To) calculated 
from Eq. (13) for (TMTSF) 2 Re0 4 at uj c t = v^O- 



The eigenvalues of the matrix (23) give the transverse 
electron dispersion relation, 



ex = 2t c cos(k z c)± ^Jcos 2 (k v b)[2t b + 4t' c cos(k z c)} 2 + E%. 

(24) 

Expanding Eq. (24) to the zeroth and first order in t c 
and comparing it with Eq. (3), we find s y ((f>) to be the 
same as in Eq. (21) and 



/(</>) = 1 ± 



t' 



2 cos 2 



t c y/cos 2 4> + {E g /2t b f 



(25) 



Only the second term in Eq. (25) generates the coeffi- 
cients A n with n^O when substituted into Eq. (12) at 
B x — 0. For Eg/tb <C 1, this term is close to a rectified 
cosine signal, as shown in Fig. 2(b), and its Fourier coef- 
ficients decay as 1/n 2 for large n. It has non-zero Fourier 
coefficients only for even n, thus a zz vs B' has peaks at 
the even Lebed angles, as shown in Fig. 5 for B' x = 
and observed experimentally in (TMTSF) 2 C10 4 4 ' 5 . Be- 
cause the second term in Eq. (25) is highly nonlin- 
ear in cos <p, it generates many harmonics and many 
Lebed peaks. However, they decay with the increase 
in n faster in (TMTSF) 2 C10 4 than in (TMTSF) 2 Re0 4 . 
Moreover, because t' c is small, the Lebed oscillations in 
(TMTSF)2C10 4 are weak, in agreement with the obser- 
vations in Rcfs. 4,5,39. As discussed in Sec. Ill, the DKC 
oscillations are controlled by the coefficient Aq(B' x ) in 
Eq. (13). The first term in Eq. (25) gives the main con- 
tribution to Aq(B' x ), proportional to Jq(B x ). Thus, the 
DKC oscillations are relatively strong in (TMTSF) 2 C10 4 , 
as observed in Refs. 8,39, although they are somewhat re- 
duced by the second term in Eq. (25). 

We conclude that the different types of anion ordering 
in (TMTSF) 2 Re0 4 and (TMTSF) 2 C10 4 can indeed ex- 
plain the characteristic features of AMRO in these ma- 
terials. In (TMTSF) 2 Re0 4 , the Lebed oscillations are 
strong and numerous, but the DKC oscillations are very 
weak. In (TMTSF) 2 C10 4 , the Lebed oscillations are nu- 



FIG. 5: (Color online) Plot of (Ji/ctq vs. B' y at B' x — 0, shown 
in the linear (left) and logarithmic (right) scales. 02 is the 
contribution to <r zz in Eq. (13) from the second term in Eq. 
(25), and a' = ao(t' c /t c ) 2 . 



merous, but weak, whereas the DKC oscillations are rel- 
atively strong. On the other hand, there is no anion 
superstructure in (TMTSF) 2 PF 6 . This material exhibits 
a few but strong Lebed oscillations and partially sup- 
pressed DKC oscillations. 



VI. INTERBAND TUNNELING IN A STRONG 
MAGNETIC FIELD PARALLEL TO THE LAYERS 

Folding of the Brillouin zone due to anion ordering 
produces two branches (or two bands) of the electron 
dispersion, which we label by the index a = ± according 
to the sign in Eq. (21). The Fermi surfaces of the two 
bands, obtained from Eq. (2), are shown by the two solid 
lines in Fig. 6 for E g /tb = 0.1. (This picture is for the 
Fermi-surface sheets near +fcp.) 

In this section, we study the interlayer conductivity 
in a strong magnetic field (B x , B y , 0) parallel to the lay- 
ers. We use the formalism developed in Refs. 35-37 and 
calculate a zz between just two layers, i.e., for a bilayer. 
Assuming that t c is very weak, one can argue that, in 
the lowest order in t c , the interlayer conductivity of a 
bulk multilayer crystal is determined by the interlayer 
conductivity between a pair of layers 56 . 

The tunneling Hamiltonian between layers 1 and 2 is 

H c = t c J $\ (r ) ^1 (r ) e ix(r) d 2 r + H.c. , (26) 

6C 

X (r) = jA z (r), A z (r) = B x y - B y x, (27) 

where -01,2 arc the electron destruction operators in lay- 
ers 1 and 2. Here A z is the vector potential of the in- 
plane magnetic field, and x(r) is the corresponding gauge 
phase accumulated in the process of tunneling across the 
interlayer spacing c. Substituting Eq. (27) into Eq. (26) 
and using momentum representation in the [x, y) plane, 




FIG. 6: (Color online) Fermi surfaces of two adjacent layers 
shifted by the vector q of Eq. (28) due to an in-plane magnetic 
field. The Fermi surfaces for each layer (the solid lines and 
the dashed lines) consist of two bands separated by the gap 
2E 9 /vf due to anion ordering and labeled + and — . 



we observe that the in-plane wave vector of the electron 
changes from k to fc + q in the process of tunneling 37 , 
where the vector q is 



(qx,q y ) = —(B y ,-B x ). 



(28) 



Thus, the Fermi surfaces of the second layer are shifted 
by the vector q relative to the Fermi surfaces of the first 
layer as shown by the two dashed lines in Fig. 6. A 
similar picture was discussed for closed Fermi surfaces in 
semiconducting bilayers in Refs. 35,57-59. 

The interlayer conductivity er"/ between the bands a 
and (3 is given by the following expression 35,60 



<*f}\ 



2 S(k,E F )S(k + q,E F ), (29) 



where M af 3 — (?pa\k + q)\ipp(k)) is the scalar product 
between the in-plane electron wave functions belonging 
to adjacent layers. These matrix elements are discussed 
in more detail in Appendix A. The total interlayer con- 
ductivity is the sum over all bands a zz =^2 a a & zz ■ The 
function S(k, E F ) is the spectral density of the in-plane 
electron Green's function evaluated at the Fermi energy 
E F as a function of the wave vector fc 35 - 60 



S(k,E F ) 



2T 



[E F -e(k)] 2 +T 2: 



(30) 



where T = H/2t is the relaxation rate, and e(k) is the 
electron dispersion within the layer. 

When r is small, i.e., when the electron quasiparticles 
have a long lifetime r, the spectral function (30) can be 
replaced by a delta function, S(k 7 E F ) w 2n5[E F — e(k)]. 
Substituting this expression into Eq. (29), we find 



«/3 _ 



Tin 



■ J J dk y dk x 5[hvpk x + a2t b e y (k y b)) 



FIG. 7: (Color online) Phase diagram of interlayer tunneling 
vs. the normalized in-plane magnetic-field components B y and 
B x . Tunneling between the same and different types of bands 
is possible in the upper left and the lower right regions of the 
diagram, correspondingly, and not possible in the intermedi- 
ate region. The thin curves show the interlayer conductivity 
a zz calculated using Eq. (34) as a function of B y for several 
values of B x for the superstructure of (TMTSF^ReOzi. 



x 5[hv F (k x - q x ) + P2t b e y (k y b - q y b)], 



(31) 



where the matrix element M a p is evaluated at the points 
where both delta functions are satisfied. Integrating 
Eq. (31) over k x , we find 



T af3 



eHlc\M aP \ 2 
nfi 2 vp 



dky 6\g a p(ky)], (32) 



where the function g a f3(k y ) is 

g a ii{k y ) = v F q x + 2t b [ae y {k y b) - f3e y {k y b - q y b)\. (33) 
Taking the integral (32), we find 



af} = eH\c v \M a p(ky)\ 2 
zz Trh 2 v F ^\dg af3 /dky\' 



(34) 



where the sum is taken over the points k y where the 
equation g a p{ky) = is satisfied. Notice that the relax- 
ation time r drops out from Eq. (34), so a zz should be 
temperature-independent in a strong parallel magnetic 
field 35 . 

Equation (31) shows that a non-zero contribution to in- 
terlayer conductivity comes from the points where both 
delta functions are satisfied, i.e., the initial and final 
states belong to the Fermi surfaces of different layers. 
Geometrically, these are the intersection points ky 1 ^ and 

7 (2) 

k y of the solid and dashed lines in Fig. 6. Depending 
on which Fermi surfaces intersect in Fig. 6, electrons can 
tunnel between different bands a, 8 = ± in the folded 
Brillouin zone. The equation g a f3(k y ) — has solutions 
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FIG. 8: (Color online) The same as in Fig. 7, but for the 
superstructure of (TMTSF) 2 C10 4 . 



only in some regions of the (q x ,q y ) space, as shown by 
the thick solid lines in Figs. 7 and 8. Above the diagonal 
line in Figs. 7 and 8, the interlayer tunneling is possible 
only between the bands of the same type a = (3. If q x 
exceeds the threshold value, 



Tt- v fQx — B y ecvp > 2E„ 



(35) 



the interlayer tunneling between different bands, a = 
— (3, becomes possible in the lower right region in Figs. 
7 and 8. No interlayer tunneling is possible in the inter- 
mediate region in Figs. 7 and 8, where the shifted Fermi 
surfaces in Fig. 6 do not cross. The boundaries of the re- 
gions are determined by the condition that the displaced 
Fermi surface touches the other one. 

The plots of the interlayer conductivity a zz , calculated 
from Eq. (34), are shown in Figs. 7 and 8 as functions 
of By oc q x for several fixed values of B x oc q y . We ob- 
serve that the interlayer conductivity vanishes in the in- 
termediate region and has peaks at the boundaries. The 
peaks originate from the increase in the phase volume in 
the integral (31) when the two Fermi surfaces touch each 
other. Figure 7 corresponds to the anion superstructure 
of (TMTSF) 2 Re0 4 . We observe that, when the magnetic 
field is applied along the y axis {B x = 0), the interlayer 
conductivity o zz (B y ) is strongly suppressed until B y ex- 
ceeds the threshold, and then a zz increases sharply. The 
value of Eg can be determined from the measured thresh- 
old field B y via Eq. (35). Figure 8 corresponds to the an- 
ion superstructure of (TMTSF)2C104. In this case, the 
eigenfunctions of different bands a = —(3 are orthogo- 
nal, so the matrix element M |_ vanishes for q y = (see 

Appendix A). Thus, in order to get a nonzero interlayer 
conductivity in (TMTSF)2C104, it is necessary to have 
a non-zero component B x ^ 0, so that q y ^ 0. 

According to the measurements in Ref. 18, the Fermi 
velocity in (TMTSF) 2 C10 4 is v F « 10 5 m/s. Substi- 
tuting this value and the interlayer distance c = 1.35 
nm (Ref. 1) into Eq. (35) and using the maximal sta- 



FIG. 9: (Color online) The in-plane Fermi surface of k- 
(ET)2Cu(NCS)2- The a and (3 branches of the Fermi surface 
are separated by the distance Afc in the momentum space. 



tionary field of 45 T available at NHMFL in Tallahassee, 
we find the maximal anion gap 2E g w 70 K that can 
be probed using this method. Various estimates of E g 
are reviewed in Ref. 61. Refs. 27,62 estimated E g as 
40 -7- 50 K, so the field of 45 T may be sufficient to exceed 
the threshold (35) at the ambient pressure. The exper- 
iment can be also performed in pulsed fields or under 
pressure, where the anion superstructure is progressively 
suppressed 63 . Measurements of the interlayer conductiv- 
ity using pulsed magnetic fields of 46 T were performed 
in (TMTSF) 2 C10 4 64 , but the field was applied close to 
the x axis, rather than to the y axis, as required for our 
effect. 

A similar analysis can be also applied to the mate- 
rial k-(ET)2Cu(NCS)2, whose in-plane Fermi surface is 
shown in Fig. 9. The separation Afc between the a and (3 
branches of the Fermi surface can be measured by apply- 
ing an in-plane magnetic field in the horizontal direction 
in Fig. 9. This field shifts the Fermi surface of one layer 
by the vector q shown in Fig. 9. The threshold magnetic 
field, at which the a branch in one layer starts to touch 
the (3 branch in the other layer, can be calculated from 
Eq. (28). Using Ak = 0.17 nm -1 and the interlayer dis- 
tance c = 1.62 nm 1,65 , we estimate that the threshold 
magnetic field is of the order of 430 T, which is beyond 
the current experimental capabilities. 



VII. CONCLUSIONS 

We have shown that the modifications of the electron 
dispersion due to the anion ordering in (TMTSF)2Re0 4 
and (TMTSF)2C10 4 generate effective tunneling am- 
plitudes between many distant chains. These ampli- 
tudes cause peaks in the interlayer conductivity cr zz 
at many Lebed magic angles (1). The different wave 
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vectors of the anion ordering, Q = (0,1/2,1/2) in 
(TMTSF) 2 Re0 4 and Q = (0, 1/2,0) in (TMTSF) 2 C10 4 , 
result in the odd and even Lebed magic angles, as ob- 
served experimentally 40,63 . Our theory also explains why 
the Lebed oscillations are strong and the DKC oscil- 
lations are weak in (TMTSF)2Rc04, and vice versa in 
(TMTSF) 2 C10 4 , as observed experimentally 39 . 

When a strong magnetic field is applied parallel to the 
layers and B y exceeds a certain threshold, then inter- 
layer tunneling between different branches of the Fermi 
surface, produced by folding of the Brillouin zone, should 
become possible. This effect would be observed as a sharp 
increase in interlayer conductivity. It can be utilized for 
a direct measurement of the anion gap E g . Theoretical 
description of this effect required a quantum-mechanical 
treatment of the wave functions confined to different lay- 
ers and cannot be achieved within the framework of qua- 
siclassical electron orbits on a warped Fermi surface. 

Experimental observation of the high number of magic 
angles (up to 21 in Ref. 40) demonstrates a very high level 
of quantum coherence achieved in the Q1D organic con- 
ductors at low temperatures. This is remarkable given 
that the (TMTSF) 2 X materials have strong electron in- 
teractions. In different parts of their rich phase diagram, 
these materials have the Mott insulating phase and other 
exotic phases 1 ' 2 . It would be very interesting to study 
what happens to AMRO when the system is driven to- 
ward the Mott state using pressure or other variables. 

We point out that the theory of the angular mag- 
netoresistance oscillations (AMRO) in Q1D conductors 
is equivalent to the mathematically description of the 
Mach-Zehnder interference in a driven superconducting 
qubit 49 ~ 53 and of laser cooling in ion traps 54 . The simi- 
larity in the behavior of these systems demonstrates that 
quantum coherence in the Q1D organic conductors at low 
temperatures is as high as in the superconducting qubits 
and ion traps, which are actively considered for applica- 
tions in quantum computing and quantum information. 
Thus, the physics of Q1D conductors may have applica- 
tions in quantum engineering well beyond the domain of 
solid-state material science. 



Acknowledgments 

V.M.Y. is grateful for the discussions with W. Kang, 
S. Hill, M.J. Naughton, S. Uji, W.D Oliver, and S. Ash- 
hab, and for the e-mail communications with A.G. Lebed 
and L. Levitov. 



APPENDIX A: CALCULATION OF THE 
MATRIX ELEMENTS 

In this appendix, we calculate the matrix elements of 
interlayer tunneling introduced in Eq. (29). 

In the case of (TMTSF) 2 C10 4 , the interlayer tunneling 
with the amplitude t c occurs between the chains of the 



same type, as shown in Fig. 1(b) 56 . The in-plane Hamil- 
tonians of two adjacent layers are given by the same ex- 
pression, 

a=(» \ h , ^ «»(*«*) y (Ai) 

^ 2t b cos(k y b) -E g J y ' 

The eigenvalues A± and the eigenvectors \ip±) of the 
Hamiltonian (Al) are 

A± = ±y/[2t b cos(k y b)] 2 +E*, (A2) 

H±{k y )) = ^(\±+ E g , 2t b cos k y b), (A3) 

N a = y/[2t b cos(k y b)} 2 + (\ a +Eg) 2 . (A4) 

The matrix elements of tunneling arc proportional to the 
scalar products of the wave functions in adjacent layers: 

M__ = ^.{k v + q y )\i,-{k v )), (A5) 
M++ = {i, + {k y + q y )\i>+{k y )) (A6) 

for tunneling between the same kinds of bands and 

M_ + = ^_(k y + q y )\i>+(k y )) (A7) 

between different kinds of bands. It is clear from Eq. (A7) 

that M y vanishes for q y = because \ip + (k y )) and 

\ij)_{kyj) are orthogonal. 

In the case of (TMTSF) 2 Rc0 4 , the inter-layer tunnel- 
ing with the amplitude t c occurs between the chains of 
different types. The in-plane Hamiltonian of one layer 
has the form (Al), whereas the sign of E g is reversed in 
the Hamiltonian H' of another layer 

H '=( 0i - E « M 2 WM)y (A8) 

\2t b cos(kyb) E g J y ' 

The eigenvalues of H 1 are the same as in Eq. (A2), but 
the corresponding eigenvectors are different, 

W±(k v )) = ^r(X±- Eg, 2t b cos k y b), (A9) 

N' + = N_, N'_=N + . (A10) 

The scalar products of the wave functions in the adjacent 
layers now are 

M__ = W_(k y + qy)\iP-(k y )), (All) 
M ++ = W + {k y + q y )W + (k y )) (A12) 

for the same kinds of bands and 

M_+ = y/_(ky + q y M+(k y )) (A13) 

for different kinds of bands. Now M |_ does not van- 
ish for q y = 0, because \if)+(k y )) and \ip'_(k y )} are not 
orthogonal. 



10 



1 T. Ishiguro, K. Yamaji, and G. Saito, Organic Supercon- 
ductors (Springer, Berlin, 1998). 

2 Physics of Organic Superconductors and Conductors, 
edited by A. G. Lebed (Springer, Berlin, 2008). 

3 A. G. Lebed, JETP Lett. 43, 174 (1986); A. G. Lebed and 
P. Bak, Phys. Rev. Lett. 63, 1315 (1989). 

4 T. Osada, A. Kawasumi, S. Kagoshima, N. Miura, and 

G. Saito, Phys. Rev. Lett. 66, 1525 (1991). 

5 M. J. Naughton, O. H. Chung, M. Chaparala, X. Bu, and 
P. Coppens, Phys. Rev. Lett. 67, 3712 (1991). 

6 W. Kang, S. T. Hannahs, and P. M. Chaikin, Phys. 
Rev. Lett. 69, 2827 (1992); E. I. Chashechkina and 
P. M. Chaikin, Phys. Rev. Lett. 80, 2181 (1998); 

D. G. Clarke, S. P. Strong, P. M. Chaikin, and 

E. I. Chashechkina, Science 279, 2071 (1998). 

7 K. Behnia, M. Ribault, and C. Lenoir, Europhys. Lett. 25, 
285 (1994). 

8 G. M. Danner, W. Kang, and P. M. Chaikin, Phys. Rev. 
Lett. 72, 3714 (1994). 

9 G. M. Danner and P. M. Chaikin, Phys. Rev. Lett. 75, 
4690 (1995). 

10 H. Yoshino, K. Saito, K. Kikuchi, H. Nishikawa, 
K. Kobayashi, and I. Ikemoto, J. Phys. Soc. Jpn. 64, 
2307 (1995); H. Yoshino, K. Murata, T. Sasaki, K. Saito, 

H. Nishikawa, K. Kikuchi, K. Kobayashi, and I. Ikemoto, 
J. Phys. Soc. Jpn. 66, 2248 (1997); H. Yoshino, K. Saito,, 
H. Nishikawa, K. Kikuchi, K. Kobayashi, and I. Ikemoto, 
J. Phys. Soc. Jpn. 66, 2410 (1997). 

11 T. Osada, S. Kagoshima, and N. Miura, Phys. Rev. Lett. 
77, 5261 (1996). 

12 A. G. Lebed and N. N. Bagmet, Phys. Rev. B 55, R8654 

(1997) . 

13 I. J. Lee and M. J. Naughton, Phys. Rev. B 57, 7423 

(1998) . 

14 I. J. Lee and M. J. Naughton, Phys. Rev. B 58, R13343 
(1998). 

15 W. Kang, T. Osada, Y. J. Jo, and H. Kang, Phys. Rev. 
Lett. 99, 017002 (2007). 

16 W. Kang, Phys. Rev. B 76, 193103 (2007). 

17 A. Ardavan, J. M. Schrama, S. J. Blundell, J. Single- 
ton, W. Hayes, M. Kurmoo, P. Day, and P. Goy, Phys. 
Rev. Lett. 81, 713 (1998); A. E. Kovalev, S. Hill, and 
J. S. Quails, Phys. Rev. B 66, 134513 (2002); Y. Oshima, 
M. Kimata, K. Kishigi, H. Ohta, K. Koyama, M. Mo- 
tokawa, H. Nishikawa, K. Kikuchi, and I. Ikemoto, Phys. 
Rev. B 68, 054526 (2003); A. E. Kovalev, S. Hill, S. Taka- 
hashi, T. N. Dhakal, S. Takasaki, J. Yamada, H. Anzai, 
and J. S. Brooks, J. Appl. Phys. 93, 8665 (2003). 

18 S. Takahashi, S. Hill, S. Takasaki, J. Yamada, and H. An- 
zai, Phys. Rev. B 72, 024540 (2005). 

19 S. Hill, Phys. Rev. B 55, 4931 (1997). 

20 S. J. Blundell and J. Singleton, Phys. Rev. B 53, 5609 
(1996); S. J. Blundell, A. Ardavan, and J. Singleton, Phys. 
Rev. B 55, R6129 (1997). 

21 S. Takahashi, A. Betancur-Rodiguez, S. Hill, S. Takasaki, 
J. Yamada, and H. Anzai, J. Low Temp. Phys. 142, 315 
(2006). 

22 S. Hill and S. Takahashi in Ref. 2, p. 457. 

23 W. Wu, P. M. Chaikin, W. Kang, J. Shinagawa, W. Yu, 
and S. E. Brown, Phys. Rev. Lett. 94, 097004 (2005). 

24 T. Osada, S. Kagoshima, and N. Miura, Phys. Rev. B 46, 



1812 (1992). 

25 A. G. Lebed and M. J. Naughton, Phys. Rev. Lett. 91, 
187003 (2003). 

26 A. G. Lebed, H. I. Ha, and M. J. Naughton, Phys. Rev. B 
71, 132504 (2005). 

27 H. I. Ha, A. G. Lebed, and M. J. Naughton, Phys. Rev. B 
73, 033107 (2006). 

28 T. Osada, N. Kami, R. Kondo, and S. Kagoshima, Synth. 
Met. 103, 2024 (1999). 

29 H. Yoshino and K. Murata, J. Phys. Soc. Jpn. 68, 3027 
(1999). 

30 H. Yoshino, A. Oda, T. Sasaki, T. Hanajiri, J. Yamada, 
S. Nakatsuji, H. Anzai, and K. Murata, J. Phys. Soc. Jpn. 
68, 3142 (1999). 

31 K. Kobayashi, M. Saito, E. Ohmichi, and T. Osada, Phys. 
Rev. Lett. 96, 126601 (2006). 

32 W. Shockley, Phys. Rev. 79, 191 (1950). 

33 R. G. Chambers, Proc. Phys. Soc, London, Sect. A 65, 
458 (1952). 

34 J. M. Ziman, Principles of Theory of Solids (Cambridge 
University Press, Cambridge, England, 1972). 

35 R. H. McKenzie and P. Moses, Phys. Rev. Lett. 81, 4492 
(1998); P. Moses and R. H. McKenzie, Phys. Rev. B 60, 
7998 (1999); U. Lundin and R. H. McKenzie, Phys. Rev. 
B 70, 235122 (2004). 

36 T. Osada and M. Kuraguchi, Synth. Met. 133-134, 75 
(2003); T. Osada, K. Kobayashi, and E. Ohmichi, Synth. 
Met. 135-136, 653 (2003); T. Osada, Physica E 12, 
272 (2002); T. Osada, M. Kuraguchi, K. Kobayashi, and 
E. Ohmichi, Synth. Met. 18, 200 (2003). 

37 B. K. Cooper and V. M. Yakovenko, Phys. Rev. Lett 96, 
037001 (2006). 

38 W. Wu, I. J. Lee, and P. M. Chaikin, Phys. Rev. Lett. 
91, 056601 (2003); N. P. Ong, W. Wu, P. M. Chaikin, and 
P. W. Anderson, Europhys. Lett. 66,579 (2004); W. Wu, 
N. P. Ong, and P. M. Chaikin, Phys. Rev. B 72, 235116 
(2005); E. S. Choi, J. S. Brooks, H. Kang, Y. J. Jo, and 
W. Kang, Phys. Rev. Lett. 95, 187001 (2005); M. S. Nam, 
A. Ardavan, W. Wu, and P. M. Chaikin, Phys. Rev. B 74, 
073105 (2006); W. Wu and P. M. Chaikin, Phys. Rev. B 
76, 153102 (2007). 

39 W. Kang, Y. J. Jo, and H. Y. Kang, J. Phys.: Conf. Ser. 

51, 355 (2006). 

40 H. Kang, Y. J. Jo, S. Uji, and W. Kang, Phys. Rev. B 68, 
132508 (2003). 

41 E. I. Chashechkina and P. M. Chaikin, Phys. Rev. B 65, 
012405 (2001). 

42 K. Maki, Phys. Rev. B 45, R5111 (1992). 

43 A. G. Lebed, N. N. Bagmet, and M. J. Naughton, Phys. 
Rev. Lett. 93, 157006 (2004). 

44 A. G. Lebed, N. N. Bagmet, and M. J. Naughton, J. Phys. 
IV (France) 114, 77 (2004). 

45 V. M. Yakovenko, Europhys. Lett. 3, 1041 (1987); Sov. 
Phys. JETP 66, 355 (1987). 

46 V. M. Yakovenko, Phys. Rev. Lett. 61, 2276 (1988). 

47 A. G. Lebed, Phys. Rev. Lett. 95, 247003 (2005). 

48 N. Joo, P. Auban-Senzier, C. R. Pasquier, S. Yonezawa, 
R. Higashinaka, Y. Maeno, S. Haddad, S. Charfi-Kaddour, 
M. Heritier, K. Bechgaard and D. Jerome, Eur. Phys. J. B 

52, 337 (2006). 

49 W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. Lev- 



11 



itov, and T. P. Orlando, Science 310, 1653 (2005). 
D. M. Berns, W. D. Oliver, S. O. Valenzuela, A. V. Shytov, 
K. K. Berggren, L. S. Levitov, and T. P. Orlando, Phys. 
Rev. Lett 97, 150502 (2006). 

M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and 
P. Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 
A. Izmalkov, S. H. W. van der Ploeg, S. N. Shevchenko, 
M. Grajcar, E. Il'ichev, U. Hubner, A. N. Omelyanchouk, 
and H.-G. Meyer, Phys. Rev. Lett. 101, 017003 (2008). 
S. Ashhab, J. R. Johansson, A. M. Zagoskin, and F. Nori, 
Phys. Rev. A 75, 063414 (2007). 

R. G. DeVoe, J. Hoffnagle, and R. G. Brewer, Phys. Rev. 
A 39, 4362 (1989); R. Bliimel, C. Kappler, W. Quint, and 
H. Walther, Phys. Rev. A 40, 808 (1989). 
Eq. (12) in Ref. 37 gives a similar expression for a zz in 
the case of multiple tunneling amplitudes ti. However, it 
contains an error and should be replaced by the correct 
Eqs. (13) and (16) given here. 

In Sees. IV and VI, we ignore the tunneling amplitude t' c 
introduced in Sec. Ill and consider only the main ampli- 
tude t c . The effects discussed in Sees. IV and VI can be 
obtained already at t' c — 0. 

J. P. Eisenstein, T. J. Gramila, L. N. Pfeiffer, and 
K. W. West, Phys. Rev. B 44, 6511 (1991). 
J. A. Simmons, S. K. Lyo, J. F. Klem, M. E. Sherwin, and 
J. R. Wendt, Phys. Rev. B 47, 15741 (1993). 



V. M. Yakovenko and B. K. Cooper, Physica E 34, 128 
(2006). 

G. D. Mahan, Many-Particle Physics, 2nd ed. (Plenum, 
New York, 1990). 

S. Haddad, S. Charfi-Kaddour, M. Heritier, and R. Ben- 
naceur, Phys. Rev. B 72, 085104 (2005). 
S. Uji, T. Terashima, H. Aoki, J. S. Brooks, M. Tokumoto, 
S. Takasaki, J. Yamada, and H. Anzai, Phys. Rev. B 53, 
14399 (1996). 

W. Kang, S. T. Hannahs, and P. M. Chaikin, Phys. Rev. 
Lett. 70, 3091 (1993); H. Shinagawa, S. Kagoshima, T. Os- 
ada, and N. Miura, Physica B 201, 490 (1994); H. Shina- 
gawa, S. Kagoshima, T. Osada, and N. Miura, Synth. Met. 
70, 759 (1995); E. I. Chashechkina and P. M. Chaikin, 
Phys. Rev. B 56, 13658 (1997). 

H. Yoshino, Z. Bayindir, J. Roy, B. Show, H.-I. Ha, 
A. G. Lebed, M. J. Naughton, K. Kikuchi, H. Nishikawa, 
and K. Murata, J. Phys.: Conf. Ser. 51, 339 (2006); 
H. Yoshino, Z. Bayindir, J. Roy, B. Shaw, H.-I. Ha, 
A. Lebed, and M. J. Naughton, J. Low Temp. Phys. 142, 
323 (2007). 

P. A. Goddard, S. J. Blundell, J. Singleton, R. D. McDon- 
ald, A. Ardavan, A. Narduzzo, J. A. Schlueter, A. M. Kini, 
and T. Sasaki, Phys. Rev. B 69, 174509 (2004). 



